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1 Introduction 

This report describes continued efforts to evaluate a breadboard of a Planetary Imaging Fourier 
Transform Spectrometer (PIFTS). The PIFTS breadboard was developed under prior PIDDP fund- 
ing That effort is described in the final report for NASA Grant NAG5-6248 and in two conference 
papers (Sromovsky et al. 2000; Revercomb et al. 2000). The PIFTS breadboard was designed for 
near-IR (1-5.2 pm) imaging of planetary targets with spectral resolving powers of several hundred 
to several thousand, using an InSb detector array providing at least 64x64 pixels imaging detaib 
The major focus of the development effort was to combine existing technologies to produce a small 
and low power design compatible with a very low mass flyable instrument. The objective of this 
erant (NAG5- 10729) was further characterization of the breadboard performance, me u mg m- 
tercomparisons with the highly accurate non-imaging Advanced Emitted Radiance Interferometer 
(AERI) (Revercomb et al. 1994; Best et al. 1997). Before describing new characterization results, 
we first provide needed background information describing the PIFTS breadboard. 


2 Description of the PIFTS breadboard 

2.1 Overview. 

The PIFTS breadboard instrument consists of (1) a miniature interferometer subsystem, (2) a 
dewar-mounted InSb detector array with readout electronics and a camera control computer (3) a 
standard IR imaging lens, (4) an optical bench that holds the subsystems in alignment, ( ) immg 
and triggering electronics, and (6) external sources for testing system performance. The functional 
relationships of these components are illustrated in Fig. 1. Shaded blocks indicate components 
procured from commercial vendors. The lightest shading indicates the camera system with darker 
shades marking the Bomem interferometer and the camera lens. The remaining elements were 
either borrowed or designed and fabricated at SSEC. The scene mirror, blackbody sources, standard 
lamps, and controllers are elements that we borrowed from other SSEC projects. 
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Figure 1: PIFTS breadboard functional block diagram 
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The interferometer package and the camera system are 
the left is the Bomem interferometer package, which is 5 in x 5 m x 3.5 in. The dewax on the rignt 
, - an t cv. detector array and has a commercial IR camera lens (black housing) attached on 
rrrrro^"anTri« ^ package attached on the tight Other support 
el^ctronic^and^alibration hardware are illustrated in Fig. 3. The blackbody calibration sources 
(A and B in the figure) are of the same type as used in the AERI calibration system ( es e a. 

1997). 



Figure 2: PIFTS modules from left to right: interferometer assembly (box), lens, detector dewa 
and readout electronics, mounted on aluminum plate. 



Figure 3- PIFTS support electronics: Heated (A) and ambient (B) blackbody cavities, blackbody 
controller (C), interferometer scan drive and metrology electronics (D). t.mer/con, roller board (E), 
and camera control computer/data system (F). 


2.2 Interferometer design. 

The details of the interferometer design are illustrated in the photo and drawing in Fig. 4. The input 
beam enters through an aperture on the left face, is divided into two beams by a CaF 2 beamsplitter 
f A at 4^° to the input axis) which are recombined after reflection from the two corner cubes, 

ind exits near the bottom left (at 90" to the input). The shaded wishbone supports the two corner 
cube reflectors that provide differential delay scan when the wishbone pivots under control 
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Fieure 4- PIFTS interferometer assembly interior mechanism photograph (left) and mechanical 
drawing top view (right), for which the input from the scene is at the left and the output to the 

lens is at the bottom. 


two voice-coil drivers. In a standard Michelson configuration, there is one fixed mirror and a secon 
movable mirror to produce an optical path difference between the two beams Our use of corner 
cubes removes the need for precise alignment of the mirrors. The optical delay 18 * 

monitoring the fringe intensity of a diode laser, using a retarder plate on one path so that detectors 
with^polarizere can keep track of absolute fringe counts in either scan direction^ The scan drive 
system uses the detected laser fringe signal as feedback to maintain a constant delay scan and 
nrovide zero crossing pulses that are used to control the camera integration timing. 

PI ° At ^tThighest * resolution, the data volume produced by PIFTS is substantial. For our laser 
wavelength of 852 nm (11,737 cm' 1 ) and two samples per fringe (one at upwaxc 
one at each downward crossing), the max optical delay scan from -0.4 cm to +0 4 cm generates 
18,779 frame synch pulses per double-sided interferogram. For 64x64 pixe s an ' J, 

per pixel we obtain 153.8 MB of image data per scan. At a scan rate of 500 laser fnnges/sec 
M000 image frames/sec), the 0.8-cm scan takes 18.79 seconds and provides a maximum unapodized 
spectraTr^olution of l/(2x0.4 cm) = 1.25 cm" >, corresponding to very respectable reso ving powers 
of 1600 at 5 pm, and 8,000 at 1 pm. For applications requiring resolving powers of a few hundr 
(mineral or frost identification for example) delay scans can be shortened dramatically 

A relatively large beam-splitter is used to accommodate the expanding FOV of the imag g 
lens Even so, this design has 1/4 the volume of the commercial instrument. There is only a sma 
amount of vignetting within the 64X64 detector region for the primary input beam. However, 
there k! ^gnificant vignetting of the second interferogram that is created by radiation entering from 
thl direcHon opposite to that of the scene radiation. If the optical system had expanded in both 
directions from the beam splitter, and the aperture of the lens were sufficient to admit that large 
a beam then the second interferogram would be formed mainly from the cold finger on which the 
detector is mounted. But in our configuration there is only a small region in which t[ ' c ,ack ™ 
interferogram actually is formed from the cold finger radiation (plus whatever reflects off the optical 
elements) The outer part of the field of view actually sees background radiation emitted from t 
lens mount and the interferometer case (more information on this topic is provided in Section . ). 
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2.3 Imaging detector array. 

The detector array and readout electronics were purchased . . « f ■£££*£%£ 

c i Rarhara Focalnlane. The camera uses an InSb array of 320x25b pixels, providing a 

subarray of 9,5% operable pixels snapshot ^ 

J- 

rZll?o?"rons n t tie multiplexer included whh 

S£": s Th 2 ! g££Z5 * 

isolation channels is relatively low crosstalk. 


2.4 Imaging optics. 

“eti AmoS standard IR camera lenses this lens a to provides “~l " T* 
and array coverage (longer focal length lenses have t too large an syslem whjle 

ZfPA, resulting in a slightly larger effective f /# than 
f/2.3). The selected lens provided a total diagonal field of vi 

subarray. 


3 New Characterization Results 

3.1 Delay Scan Servo System Characterization and Analysis 

. . .. r , i RIFT'S delav scan servo to external vibrations motivated us to better 

=5SsSSSrS=?S=^-S 

out an approximately circular Lissajous pattern with i the an 
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oscillator with the following equation of motion: 

{J + D 2 + CD + K m )d 0 {t) = ( K + KC 2 D + KKiD~ l )[9\(t - t d ) - 0 o (t - t d )], (1) 

where D = d/dt , D 2 = d 2 /dt 2 , D~ l = f dt, J is the moment of inertia of the wishbone, C is the 
natural damping time (with the servo turned off), K m is the mechanical spring constant, K is the 
overall servo gain (torque/unit position error), KC 2 is the servo gain constant for speed errors, 
KKi is the servo gain for integrated position error, is the input (desired) angular position, d Q 
is the output (actual) angular position, and t d is the time delay between the applied torque (the 
RH side of the equation) and the error used to compute the applied torque (9\ - 0 o ). The Laplace 
transform of this equation, ignoring boundary conditions, is given by 

( J S 2 + Cs + Km)e 0 {t) = (K + KC 2 s + KKi/s) exp (-t d s) [<?i(t) - 6 0 (t)} (2) 


where s is the complex angular frequency. Prom this we obtain the open and closed loop transfer 
functions ^ (if + KC 2 s + KKjs) exp j-t^s) ^ 


d>c = 


$0 


Js 2 + Cs + Km 
(Ks + KC 2 s? + KKi ) exp(-trfs) 


(4) 


1 + $0 Js 3 + (C + KC 2 exp (—t d s))s 2 + ( K m + K exp(— t d s))s + KK\ exp( t d s) 

There are three poles of the closed loop transfer function. The two off axis poles in the complex 
s-plane correspond to resonant frequencies, with frequency proportional to the imaginary (y) co- 
ordinate and damping time inversely proportional to the negative real (x) coordinate. The poles 
must be in the negative half complex plane of s for the control system to be stable. A positive real 
coordinate implies exponential growth instead of damping. Increasing gain (K) or the time delay 


(t d ) moves the poles towards instability. 

We carried out a number of tests to establish mechanical properties of the drive mechanism, 
including natural spring constant, and natural resonant frequencies. We also measured parameters 
of the servo control system, including control constants for position, derivative, and integral feed- 
back terms. We found that the most important parameter limiting servo gain was the time delay 
between a servo error and the application of a correction signal. About 1 ms of delay occurred in 
the DSP used to compute the error signal, and about 0.25 ms resulted from the inductance of the 
drive coils. The servo model that roughly agrees with the PIFTS characteristics is shown in the 
left panel of Fig. 5. The measured fractional control error is an approximation based on measuring 
the maximum amplitudes of perturbations that did not result in loss of servo control. For these 
measurements the interferometer assembly was placed on a foam pad and excited at a succession of 
fixed frequencies by using a solenoid to push on a corner of the housing, resulting in a combination 
of twisting and lateral acceleration perturbations, which were measured by a linear accelerome- 
ter. The accelerations were converted to displacement amplitudes by dividing by (27 rf) . We then 
normalized by the displacement at middle frequencies beyond the range of servo control. The dis- 
agreement in the location of the resonant dips of the model and the measurements is probably due 
to the fact that the foam pad also has a resonance. The discrepancy at very low frequencies might 
also be due to a change in the fraction of the perturbation that results in a rotational displacement. 

A model with a weaker spring and shorter delay time is shown in the right panel of that Fig. 
5. We found that the servo gain could not be increased beyond the value chosen by BOMEM 
without making the servo extremely unstable, and we in fact decided to reduce the gain somewhat 
to make it more stable. We found that servo control was limited to frequencies below about 100 
Hz, and that the best solution to improving field performance was to limit high frequency inputs 
with compliant shock mounts and wind shielding. 
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Figure 5- Servo model characteristics approximating nominal PIFTS configuration (left) and 
nrovements possible with reduced delay time, reduced mechanical spring constant, and increase 
servo gain The top plot displays the three poles in the closed loop transfer functions, the second 
provides open and closed loop transfer functions, the third provides the corresponding phases, and 
the bottom curve provides the control error. Approximate measured control errors for the nominal 
PIFTS configuration are shown as plotted points. 


The main reason for the high power consumption of the servo was the use of very stiff return 
springs on the wishbone assembly. The position term of the servo also acts like a mechanical 
spring except that the mechanical spring has no time delay. Thus, if we eliminate or greatly 
reduce the mechanical spring constant, we may not be able to compensate for it by increasing 
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- ** initialization Spence used to establish servo lock 
a ‘ W^rlXTonancez in the wishbone structure, although at relatively high frequencies that 

seems to be associated with a flexing of the corner cube mount. 

3.2 Delay scale calibration for each image pixel. 

i 7 ■ in the PIFTS interferogram data cube has a slightly different delay scale associated 

lot : - "dtrir" P tL a ;iii"he 

:ni A '"wirs 

SUm X 2 = E {L(v)aeri - L{v)piFTs) 2 /a(v) 2 ’ ^ 

where L denotes radiance spectrum, and o’ is the expected — ^ 

We expect the radiance difference to be dominated 

r eS-d PIFTS radiance errors , 

biTless than Sd ~ ^ 

rwbiirs sr =£ :“f «==£ 

^rce1ht"„ral" unit ™ taken as the uncertainty in ^ For each of 1 25 
on the focal plane array the fractional uncertainty in A m0I was within the range • > 

0 0 Tfte°r“ 'Lm^b^ best-lit values of A e „ we then used non-linear regression to constrain 
the va“r model of A«„ . We can express A. f as a function of pixel pos.t.on (,. y) in terms 

of four adjustable parameters: 

Aeff(z, y) = X max cos(arctan([(x - x Q f + {y ~ yo) 2 ] 1/2 //p))> 

, / \ are the opt i c ax is coordinates in pixels, A max is the laser wavelength at the optic 

Li"' lit U the focal length in pixel units. The parameters of the model were individually varied 
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Figure 6: y* as a lunction oi a scaie i^iui «.uc» - — “ QO >, 

laser wavelength of 852 nm, shown for average centered at image coordinates (64, 32) 


as a function of a scale factor that is the ratio of trial laser wavelength to the nominal 

. r. i imrtrrn r'nnrrlinat.Pfi 1 . 


step by step over a range of values. At each value, the other parameters were adjusted to obtain 
the smallest value of * 2 at that value of the parameter of interest, 

rescaled so that the minima were equal to the number of degrees of freedom (21 in this case) and 
the value change producing a * 2 increase of 1 defined the uncertainty for that parameter. The 

results were as follows: 


x 0 = 67.084 ±0.029 pixels 
y 0 = 42.264 ±0.045 pixels 
Amai / 852 nm = 1.0007797 ±1.6E-6 
(A mox = 852.6643 nm ±0.0014 nm ) 

(1/Amax = 11,727.945 cm' 1 ±0.019 cm" 1 ) 

/ = 1/(0.00120088 ±1.53e-6) = 832.72267 ±1.06 pixels 


For a pixel pitch of 30 pm, this implies a focal length of 24.98 mm ±0.03 mm, compared to the 

" OI Thus V welre a^to find the optic axis within 0.03-0.05 pixels, the on-axis wavenumber scale 
to within 0.0014 nm (1.6 parts per million), and the imaging focal length to within , T. .2 parte pm: 
thousand (which defines the angular scale of the image to the same precision). Slightly worse 
performance is expected from comparisons between PIFTS and calculated radiances. 

With fits to more points on the focal plane it should be possible to define any important ‘T®' 
distortions that may be present as well. (Most optical systems have some of distortmm 

Examples include the WFPC2 camera system on the Hubble Space Telescope and the 
imager on the NASA Infrared Telescope on Mauna Kea, which both need some correction 
distortion.) Distortion corrections may be important in mosaicing of images taken at different view 

angles. 
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3.3 Interpolation to a fixed wavenumber scale. 

It is not sufficient to establish the correct wavenumber scale for each image pixel. To facilitate 
data analysis algorithms that combine spectra from different pixels, it is best to put all spectra on 
the same wavenumber scale at the outset. Interpolation is needed, either in the spectral domain 
or in the delay domain. Resampling in the spectral domain can be done by first zero-filling 
interferograms to several times their original length (with the optimum factor depending on the 
degree of oversampling in the original interferogram), then transforming to the spectral domain 
thich produces a highly oversampled spectrum that can be easily interpolated at t he : desired 
spectral sampling. To facilitate analyses of data from different cubes, or different days for which 
the reference laser wavelength might be slightly different, it is best to put all ^ r ^ tlo “ J °" 
standard wavelength scale. The interpolation quality improves with the length of the zero fill, but 

the time to do the FFT increases as JVlog(iV). . . . 

Current PIFTS interferograms are significantly oversampled relative to the Nyqmst criterion. 

with 23,473 samples/cm, we obtain 11.7 samples per signal fringe at 5 pm and 7 samples per fringe 
at 3 pm. Thus, it is possible to do a reasonably accurate interpolation in the delay domain. We 
compared transforms of linear and cubic spline interpolations with the transform of a factor 
zero-fill interferogram, which we use as a standard. The errors of a cubic spline interpolation (Fig. 
7)Tre less than 0.1% from 1980 cm' 1 to 3100 cn>-\ and 0.02% or less m the lower noise regions 
of the spectrum from 2000 cm' 1 to 2900 cm" 1 . Linear interpolation errors are typically a factor 
of 100 worse than this. The cubic spline errors may be even smaller than our analysis as so 
far indicated In high-noise regions, the inherent smoothness of the spectrally band-limited signal 
["ontminated by random noise variations that are only band-limited by the digital averaging 
window, which is harder to reproduce with the same percentage accuracy. The overall performance 
of the cubic spline interpolation is more than adequate. 


HBB, 19N0V, 3PM, Cubic Spline Interp. Error 

— ■ * • T 1 ' 1 ’ ' ' 1 ' ' ' F ’ F 



Wavenumber 


Figure 7: Cubic spline interpolation errors, using a noisy PIFTS spectrum with the transform of 
an 8 x zero-fill padded interferogram as the spectral reference. 
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3.4 Characterization of Background Interferograms. 

Half of the radiation entering the input aperture of the interferometer is reflected back toward the 
source It is also true that half of the radiation entering from the exit aperture (and surround! g ) 
is reflected back towards the detector. This background interferogram has been modeled and 
demonstrated to have a modulation that is inverted relative to the modulation of the iterferogram 
generated by the incoming beam. A comparison of beam combinations leading to > the 
view interferogram and the background interferogram is shown in Fig. 8. Actually, for the PIF 
and AERI instruments the interferogram modulation produced by the incoming beam is nega iv 
whilethe modulation of the background interferogram is positive. This can be inferred from the 
analysis on pages 4 and 5 of Brault (1985), with one modification. For one of our recombined bea 
the beam-splitter reflection is internal, while for the other it is external, resulting in a relative 180 
phase shift between them that produces a negative modulation. In Brault “^^2^ * 
beams each have one external reflection, leading to positive modulation. Besides the differen g 
of modulation of the interfering beams entering from the rear, the background amplitude also has 
a different value, although not significantly different. Assuming that the beam-splitter reflectivi y 
R t ^lV a nlUon Ts satisfy Rb = 1- ft, we can write the relative interferogram amplitudes 

of these two sources as, 

(7) 

(8) 


I x (x) = KRtuRbTb^ -cos(2txi/ + ^i)]I»iH 
I 2 (x) = KR m R B T B + cos ( 2nxu + <k)]T 2 W 


2RbTb 

where Rm is the reflectivity ot either beam-return mirror and K is a common proportionality 
constant It can also be shown that they have the same compensation, so that they shou lave 
the same uncompensated phase errors, although the equation is written to include the possibility 

° f SO T™ ^“s^by^det^r is the sum of the two interferograms. The observed mod- 
ulated interferogram is thus the sum of two components and the transform of the interferogram 
becomes ,.00 

I(x)exp{2mvx)dx = G\(y) - G2 \v) 1 > 


/ GO 
-GO 


where we use O, to denote the signal spectrum of the target, while G, denotes the spectrum of the 
instrument. If one attempts to determine a spectrum Gfa) = |G(i/)| exp [ifa - G 1 fa) 2 ( ), 
resulting phase can be wildly variable at wavenumbers for which the two component spectra ha 
similar values, and thus lead to sign flips. At the center of the instrument 
arises mainly from the cold regions inside the dewar, and is relatively low. But off ax s pixels 
larger and llrger background contributions, which can even exceed the amplitude of the primary 
interferogram. The interference of these two interferograms leads to large phase 

instrument field of view. Complex calibration equations (Revercomb et i 1988 Sromo^kyMTO) 
are successful at removing the background contributions, and allow us to define the characteristics 
of the background as well as the responsivity of the system. The background radiance image 

be expressed as 


L 2 = 


G cLh LcGh —<t>2) 


( 10 ) 


Gk — G 

where fa is the phase of “normal” spectrum and fa is the phase of the background spectrum, and 
subscripts c and h refer to cold and hot reference observations. The real portion of the badtgroun 
radiance spectrum, converted to brightness temperature, is shown in Figure 9. The g 
average over the 2100-2200 cm' 1 spectral range. 
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beamsplitter 



Figure 8: Beam combinations leading to interferograms of the scene (left) and tokgr 
(right) Note that one of the background beams is externally reflected by the beam-splitte 
tilaver reflective coating (dark line facing the air gap), while for the primary mterferogram one 
Tin internal and one an external reflection from that coating. This is the source of the 

relative phase change between the two modulations. 



Fieure 9- Background radiance field (real part), shown as brightness temperature averaged over 
2 M 2200 7^ The dot-dash circle indicates the 6.22" unvignetted PIFTS FOV through the 
entrance^aperture. The reflected view through the exit aperture is vignetted by the background 
structure at angles determined mainly by the relatively small exit aperture. 


If there were no vignetting, and no reflection from optical elements between he interferometer 
and focal plane array, we would expect a relatively uniform brightness temperature around 77 , 

the temperature of the LN 2 -cooled array. Instead we see a retrieved brightness temperature of ^SO 
K at the center of the image, which can be obtained by a linear combination of 70% emission 
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K and 30% emission at the ambient temperature of 283 K. 


3.5 Defining data processing flow for producing calibrated spectral cubes. 

This summarizes the flow of information in our processing of PIFTS dat from raw data acquisition 
through creation of calibrated scene spectra. 

1 Measure a set of detector signals at delay steps Ax, as the interferometer is scanned from 
nominal delay -X x to +X U in which delay is measured relative to a reference laser that 
may be sampling an angular position different from that sampled by the scene interferogram 
detector. Let this set of readings be V k ,k = 0, N x - 1 at nominal delay positions x hk ,k - 

2. Re-scale the delay values to account for the off-axis angle of the detector relative to the laser 
reference, and for any other related scale changes. 


If this is a hot blackbody measurement find the zero delay position by locating the extreme 
value of the interferogram (this could be taken as the extreme value, or use a fit to the 
extreme value and its neighbors to find an interpolated peak). If this is not a hot blackbody 
measurement, then use the last determination of the zero delay position (which must have 
been done since the interferometer last lost servo lock). 


4 Interpolate all interferogram values to a standard set of delay positions so that both in er- 
ferograms and spectra for different detector elements in an image can be directly compared 
or averaged This will result in a new set of N interferogram values V(x k ) at delay values 
Xki k = o, N - 1 from -L to +L, where N < N\ and L < X\. Shifting and symmetrizing the 

delay range necessitates some loss of data. 

5 Subtract the mean and linear ramp, then flip the sign to get a modulation signal that is 
positive at zero delay. The result will be called V,(x fc ) which is some constant K times the 
modulated part of the interferogram radiance (defined in Eq. 3 of Sromovsky 2002). 


(V(x) for x = 0, . . . , L 

6. Create a re-ordered array V a (at) = j ^ for x = L, . . . , 2L 


7. Multiply by ^ to get h k — ^V^(x*), A: — 0, . . . , N 1 

8. Carry out discrete FFT to obtain H n , n = 0, . . . , N — 1 

'H n n = 0, . . . , N/2 — 1 

9. Reorder H n to get KH(u n ) = < 

where v n = 

function obtained from the convolution 

H(u) = G{ v) * 2Lsinc(2Li/), 


- - n - > > • 

H n+N n = —N/2 , . . . , -1 , where H is the complex spectral 


where 


G{ v) 


\B(v) exp(i<£„) for v > 0 

^B(-v)exp(-i(j)-v) for u < 0 


and where K is the conversion between spectral radiance units and measured units (DN or 
voltage). If this is a hot or cold blackbody measurement, save this result as the complex 
spectrum KH h or KH C . Otherwise, proceed to the next step. 
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10. If this is a complex scene spectrum, then compute a calibrated spectrum as 

L = L c + 


KH - KH ‘-(L h -L c ). 


KH h -KH c 

Set L equal to Re(L) to insure that imaginary components introduced by errors are eliminated. 
Note that this spectrum will be the true spectrum convolved with a sine function, provided 
that both and r,{v) are sufficiently slowly varying that the following approximation is valid: 

/ OO f 

^L(i/)sinc(2(i/ - v')L)dv . 

J — OO 

11. Compute responsivity using 


R = 


| KH h - KH C \ 

r, _ T. 


An internal report (Fry 2002) illustrates these steps in great detail. Further details concerning 
the^al(fc*ation process arj discussed in a paper submitted to Applied Op.ies (Sromovslty 2002). 

3.6 Noise in Raw and Calibrated Complex Spectra. 

Complex calibration is required to properly eliminate the effects of phase errors, complex back- 
ground signals, and anomalous phase errors due to beam-splitter emissions. The basis for complex 
calibration, and a discussion of the effects of complex calibration on noise is provided in a paper 
submitted to Applied Optics on 29 March 2002 (“Radiometric Errors in Complex Fourier Spec- 
trometry” by L. A. Sromovsky). The error in the calibrated radiances depends on the error in 
the calibration measurements, the uncertainty in the source radiances and the error in the scene 
radiance measurement. The relationship between these quantities is as follows. 

e\ s = mV5 2 +[ ai c + m 2 4c]((^-^)/(^-^)) 2 + ^ +m2 ^" 1((Ls_Lc)/(ih_ic))2 (U) 

where a VS , *vc, o VH , are the standard errors in measuring the detector voltage (or DN) signals 
in response to scene, cold blackbody, and hot blackbody references where m = (^-^)/(^h V e ) 
is the^ slope of the calibration line used to convert voltage (or DN) to radiance, and wher lc 
and o LH are the uncertainties in the radiance values of the two reference sources. The expected 
noise varies with scene radiance. The minimum variance in calibrated radiances is found Jorscene 
radiances between the two reference radiances; the variance grows quadratically away from the 

"importance of complex calibration is shown to be especially great when background inter 
ferogram amplitudes are comparable to primary interferogram amplitudes. Because the sign of the 
background modulation is opposite to that of the scene interferogram, phases can be wildly va Tying 
where net amplitudes are small due to cancellation, leading to poorly behaved magnitude spectra. 
Complex calibration eliminates this problem. 


3.7 Intercomparison with AERI observations. 

Preliminary comparisons were made between PIFTS and AERI observations on 19 November 2001. 
Both instruments were placed in close proximity outside our building (Fig. 1°) and arranged [to 
put the AERI field of view within the PIFTS field of view, as shown in Fig. 11. We also launched 
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a radiosonde from the top of our building to provide data needed for calculating the radiance 
field under observation. November 19 was a day with apparently clear and photometric skies, but 
generally high winds required us to use a wind shield (Fig. 10) to keep the PIFTS delay scan 
servo from occasionally losing lock. The high winds also resulted in rapidly increasing radiosonde 
range after release, and resulting loss of contact with the radiosonde at a relatively low altitude. 
Although this radiosonde only returned data up to 10.95 km, the radiance contributions above this 
level have been shown to be negligible for a clear atmosphere, thus providing sufficient atmospheric 
information to permit accurate line-by-line spectral radiance calculations. Prior observations on 15- 
16 November were less wind-perturbed, had better vertical coverage by radiosondes, but contained 
a small degree of cirrus cloud contamination, which might complicate comparisons with radiation 
transfer calculations. These observations were also noisier than those on 19 November. 



Figure 10: PIFTS and AERI configured for intercomparison test in the yard next to the Space 
Science and Engineering Center on the UW Madison campus. 


The PIFTS spectra that were used in the initial 19 November comparison were really only 
in preliminary form. Interferogram slopes, had not been removed prior to calculating spectra, 
although sample calculations indicate that the effects of the slope are Probably not significant. 
For the preliminary analysis, the wavenumber scale was obtained from fitting PIFTS observations 
averaged over the central 9 x 9 subimage to AERI observations. These results are shown in Fig. 

12 ’ PIFTS and AERI spectra from 19 November 2001 agree well with each other except between 
2300 and 2400 cm” 1 , where AERI results have greater errors due to longer path lengths to black- 
bodies in this region of extreme C0 2 opacity, and between 2080 and 2180 cm" 1 , for which there 
is no obvious explanation. Below 2100 cm- 1 both PIFTS and AERI disagree with calculations, 
perhaps because water amounts during PIFTS and AERI observations were different from those 

reported by the radiosonde. 
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Figure 11: Comparison of PIFTS and AERI fields of view during intercomparison tests on 19 
November 2001. The structure at the right is the north side of the Atmospheric, Oceanic, and 
Space Sciences building. 


Radiance calculations made use of the LBLRTM code developed by Clough et al. (1991). These 
calculations require that radiosonde data files be sanitized to avoid negative, zero, or very small pos- 
itive altitude increments that round to zero. In addition, obtaining good agreement of calculations 
with AERI and PIFTS measurements required adding a IK offset to the radiosonde temperature 
data. Although that correction works well, there may be some other way to achieve agreement 
without such a large “correction” to the radiosonde data - larger than local expertise would expect. 
However, it is clear that just modifying the near-surface structure (i.e. below the first radiosone 
observation) cannot produce agreement over a broad spectral range. . 

After more extensive analysis, following more accurate data processing procedures for putting a 
observations on the same wavenumber scale (described in Section 15 and more fully un the 23 pn 
2002 Fry report), the above intercomparison results for 19 November 2001 remain valid. However, 
analysis 16 November 2001 observations revealed that the laser wavenumber had drifted 1.1 parts 
in 10 4 With that correction we were able to obtain a PIFTS-AERI difference spectrum that 
had a mean value of only 0.004 mW/(m 2 -sr-cm- 1 ) and an RMS value of 0.018 mW/(m 2 -sr-cm- ) 
over wavenumbers from 2000 cm' 1 to 2600 cm" 1 . The PIFTS spectrum used in this difference 
was an average over a 35-pixel by 35-pixel box, taken to reduce random noise contributions. The 
systematic errors in this difference spectrum are seen to be extremely small. The reason for the 
laser wavenumber change remains to be determined. The laser diode emits at a wavenumber that 
depends on temperature. There is a servo that maintains the laser over a narrow temperature 
range although it is not clear at this point whether the servo is operating properly. 

The intercomparison test was partly supported by a grant from the University of Wisconsin 

Graduate School. 
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Figure 12: Comparison of PIFTS radiance from 19 November 2001 with simultaneous AERI cal- 
culations and LBLRTM calculations. 


3.8 Noise Measurements and Models. 

Random radiometric errors appear to be 0.02 to 0.04 mW/K-sr-cnr 1 ) for PIFTS relative to 
calculations, and for PIFTS relative to AERI, using 9-pixel by 9-pixel averages for PIFTS. Single- 
pixel noise levels would be presumably 9 times larger, or 0.18 to 0.36 mW/(m -sr-cm ), quite a 
large noise level in a region where typical atmospheric emitted radiance values are of the order of 
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1 mW/(m 2 -sr-cm l ) or less. ...... e . ,. 

The noise levels derived for single-pixel radiances, derived from the standard deviation of indi- 
vidual calibrated spectra for individual pixels relative to the mean value (as described in the final 
report on the original PIDDP grant), are approximately 0.075-0.08 mW/m -sr-cm , measured un- 
der laboratory conditions with a delay scan of ±0.174 cm" 1 (unapodized resolution of 2.87 cm). 
Because NESR is proportional to L/y/N when dwell time per sample is fixed, for a delay scan twice 
as long (unapodized resolution of 1.435 cm' 1 ), the noise is increased by y/2 to 0.106-0 113 ru This 
suggests that random detector noise contributions to a 9 x 9 average should be about .012-.013 ra- 
diance units (about 1/2 to 1/3 of the PIFTS-AERI difference observed during the intercompanson 
on 19 November), and 35 x 35 pixel averages should have negligible contributions from detector 

no ise. . • 

The observed rms PIFTS-AERI and PIFTS-LBLRTM differences are 0.02-0.04 radiance units, 

where LBLRTM are calculated radiance values using radiosonde profiles as inputs. It is hard to 

attribute much of this to AERI because the AERI noise level is supposed to be only 0.005 radiance 

units (R. Knuteson, personal communication, May 2002). 

A reasonable model of rms deviations between PIFTS and AERI is as follows: 


_2 

a tot 


°\eRI + a SYS + VpiFTs/N, 


( 12 ) 


where oaeri is the rms noise in the AERI spectra, a pifts is the rms noise in the single-pixel PIF 
spectra, N is the number of pixels averaged, and a S YS is the systematic rms error associate wit 
delay sampling errors, laser wavelength drifts, and other errors that affect spectra for all pixels in 
a given cube in the same way and are thus not reduced by averaging over many pixels. 

If we estimate AERI noise from the standard deviation of 5 successive AERI spectra relative 
to their mean spectrum, we find a value generally closer to 0.01 radiance units, rather than t le 

0 005 units suggested by Knuteson. This probably includes some variability in the atmosphere 
as well as instrumental noise, however. If we use 0.01 units for a aeri, then the remaining offset 
variance (the value at N=infinity) must be attributed to <r S ys- From the 14:55 cube on 15 NOV 
we find good agreement with the model, using a S ys = 0.026 ru, and a PIFT s= =0-243 ru the latter 
being established by the variation of o\ OT with N . We varied N by varying the size of the PIFTS 
subarray that is averaged. For this cube there was a clear linear dependence of total variance on 

1 iff For the 1315 cube on 19 November, which has much less noise, the dependence of a TOT on 
1/N does not closely follow the linear model. Rough agreement is obtained with a S ys = 0.023 ru 
a aeri = 0.01 ru, and a PIFT s= 0.16 ru. Even this lower noise level is much larger than expected. 

It also seems to be larger than estimated by other means. 

If we average all PIFTS spectra in a subarray and compute the standard deviation of individual 
pixel spectra from that mean, we find larger values for intercomparison tests than expected from 
laboratory noise measurements when adjusted for resolution differences. The closest agreement is 
for the 19 Nov 02 13:15 cube, from which we obtain a minimum noise of about 0.13 ru, compared 
to the value of 0.11 ru inferred from laboratory results discussed previously. One of the worst 
disagreements is for the 16 November 15:45 cube, from which we obtain a minimum no.se of 0.16 
ru. The variation is even larger outside the 2200-2300 cm' 1 region of the minima. For example at 
2450 cm -1 , which we expect to be close to 0.11 ru from lab results, ranges from 0.145 ru (19 Nov 
13T5) to 0.18 ru (16 Nov 1545). These standard deviations were computed for square subarrays, 

generally 35 pixels on a side. 

Even more troubling is the presence of asymmetries in variability between row and column 
directions (Fig. 13). If we use only a row of pixels to compute a mean and standard deviation, we 
get about the same results as using a square array. But if we use a column of pixels, then we n 
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Figure 13: Variation of standard deviations with wavelength and with direction of sampling. The 
spectral shape in these figures is mainly controlled by responsivity variations. Noise is much reduced 
on 19 November, but both 19 and 16 November data sets exhibit variability with column value. 


a large variation in results depending on column number. In the 2200-2300 cm -1 region we find 
values ranging from a low of about 0.1 ru to a high of 0.18 ru for the 19 November 15:00 cube, and 
from a low of about 0.125 ru to a high of 0.2 ru for the 15 November 14:55 cube. The variation 
with column seems to contain a roughly periodic component with a 4-pixel period. 

Thus it appears that the PIFTS noise level during the intercomparison tests was substantially 
higher than during earlier lab tests, and that there are substantial variations in noise level that 
have spatial patterns suggesting that the readout electronics may itself be excessively noisy. We 
have so far been unable to determine the cause of the increased and patterned noise levels obtained 
in recent testing. 

3.9 Publications. 

Some of the results described in previous sections became part of a paper submitted to Applied Optics 
in March 2002 (Sromovsky 2002). The paper is still under review. 
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